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All stationary solutions to the one-dimensional nonlinear 
Schrodinger equation under box and periodic boundary con- 
ditions are presented in analytic form. We consider the case 
of repulsive nonlinearity; in a companion paper we treat the 
attractive case. Our solutions take the form of stationary 
trains of dark or grey density-notch solitons. Real station- 
ary states are in one-to-one correspondence with those of the 
linear Schrodinger equation. Complex stationary states are 
uniquely nonlinear, nodeless, and symmetry-breaking. Our 
solutions apply to many physical contexts, including the Bose- 
Einstein condensate and optical pulses in fibers. 



I. INTRODUCTION 

The nonlinear Schrodinger equation (NLSE) models 
many phenomena observed in the recently created dilute 
gas Bose-Einstein condensates (BECs)El. In this conte: 
it is also referred to as the Gross-Pitaevski equationB^ 
The NLSE is ubiquitous. It describes wave propaga- 
tion phenomena in many systems besides the BEC, in- 
cluding optical pulses in fibergj, helical excitations of 
a vootex lineO, Bose-condensed photonaZl, and magnetic 
filmsB. It is one of a few basic equations upon which the 
modern Jjheory of integrable nonlinear systems has been 
foundedBEl. 

Many applications of the NLSE to BECs have dealt 
with ground-state properties, but there is growing inter- 
est in the possibility of generating topological excitations 
of a condensate, which may well be described by excited- 
state solutions of the NLSE. In this article we investigate 
such solutions, in the case of a one-dimensional NLSE 
for repulsive nonlinearity, subject either to box or peri- 
odic boundary conditions on a finite interval. This corre- 
sponds to a BEC with repulsive atomic pair interactions, 
which is the case that has received the most experimental 
interest; a companion paperEHI treats the case of attrac- 
tive interactions. The stationary excited states that we 
study here are related to the well-known soliton solutions 
of the time-dependent NLSE, and when perturbed give 
rise to soliton-type motionsU. Recent experiments show 
that such motions can be induced . in . B ECs by optical 
engineering of the condensate phasetHlEJ. 



Box and periodic boundary conditions are as ubiqui- 
tous as the NLSE and give physical iiisight into the so- 
lutions to more complicated potentialaiJ't£l. They model 
the potentially quasi-one-dimensional regime of a iium- 
ber of present experiments: the atom waveguidellSEZl; the 
prolate, haxrnonic trap in which is formed a cigaij-shaped 
BECIISIIj'Ej; the newly-developed toroidal trapL^I; and fi- 
nally an oblate harmonic trap with a barrier formed in 
the middle either by |a-| second spin state of the same 
atomOH or by a laseJlj, in which is formed a pancake- 
shaped BEC with the center removed. Periodic boundary 
conditions provide a first model for toroidal geometries; 
box boundary conditions are a good starting model for 
cigar-shaped geometries. 

The full spectrum of soliton solutions to the NLSE 
on the,-4iftfinite line was discovered by Zakharov and 
Shabatcj'O. These authors used the Inverse Scattering 
Transform, a method to which a great deal of mathemat- 
ical physics literature has been devotedO. To solve the 
NLSE under these new boundary conditions we have cho- 
sen instead to use straightforward, elementary methods, 
accessible to a broad spectrum of physicists and simpler 
than the Inverse Scattering Transform. 

There have been many recent applieaiipns of Zakhaxav; 
and Shabat's solutions to the-.^CElEl. ToroidalBEZi 
and cylindrical or box-shapedcSEa geometries have been 
considered. Some a|(ithopS| have solved the parabolic 
potential numericallyl23'M'EJ. However, to the best of 
our knowledge no one to-date has explored, analytically 
or otherwise, the full spectrum of bounded, stationary, 
multiple-soliton solutions to the NLSE under periodic 
and box boundary conditions. 



II. THE QUASI-ONE-DIMENSIONAL NLSE 

The nonlinear Schrodinger equation which describes 
the BEC for T <C Tc is a three-dimensional mean field 
theory. We define the BEC to be in the quasi-one- 
dimensional regime when its transverse dimensions are 
on the order of its healing length and its longitudinal di- 
mension is much longer than its transverse ones. In this 
case the ID limit of the 3D NLSELis appropriate, rather 
than a true ID mean field theoryo, as would be the case 
for a transverse dimension on the order of the atomic 
interaction length. 
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The quasi- ID limit of the 3D NLSE is impHcitly used in 
manyp^laces in the literature, a recent example being Liu 
et alSB Discussions which take into account transverse 
excitatioiis-||and various other geometries may be found 
elsewherccZla. We present a brief derivation in which we 
require that the wavefunction be approximately separa- 
ble and that its transverse parts be in the ground state. 
By requiring the transverse dimensions to be on the order 
of the healing length we ensure both separability of the 
wavefunction and transverse stability of the sokutians, as 
we have numerically illustrated in other worksr'lFi. 

We begin with the NLSE that describes a BEC of N 
atoms of mass M, confined in an external potential V{r): 



2M 



W^+g\i^{r,t)f +V{r)]^{r,t) 
= ihdt il^if, t) 



(1) 



where \ip{f,t)\'^ is the single particle density such that 
p{r,t) — N I i}{r,t) p, the coupling constant g = 
ATTfL^aN/M, and a is the s-wave scattering length for 
binary collisions between atoms. The case of repulsive 
interactions considered here corresponds to a > 0. 

y(r) is defined to be a three-dimensional rectangu- 
lar box of length L and small transverse area At- In 
the transverse directions the wavefunction is required to 
vanish on the surface of the container; in the longitudinal 
direction we require either box or periodic boundary con- 
ditions. This models the quasi-one-dimensional regime of 
many BEC experijneiits, as was mentioned in Sec. |, as 
well as jing lasersE3EZl, helical excitations of a vortex line 
or ringS, and many other physical systems for which the 
ID NLSE is a good model. 

The characteristic length scale over which the conden- 
sate density resumes its average value away from a sharp 
defect or from a perfectly confining wall is the healing 
length ^: 



e^(8^p|a|)-i/2 



(2) 



p = N/{LAf) is the mean particle density, where L is 
the longitudinal length of the confining potential and Ly 
and Lz are the transverse lengths. At = LyL^ is the 
transverse area. 

The BEC is in the quasi-lD regime when Ly and Lz 
satisfy the following criteria: Ly,Lz ^ ^ and Ly,Lz ^ 
L. The former ensures that the condensate remains 
in the ground state in the two transverse dimensions 
while the latter ensures that longitudinal excitations are 
much lower in energy than possible transverse excita- 
tions. Under these conditions one may make an adia- 
batic separation of longitudinal and transverse variables: 
V'(f,t) = (i^t)-i/2/(a;)/i(2/,z)e-'^*/^, where fix) and 
z) are dimensionless and the time dependence of a 
stationary state has been assumed, /i being the chemical 
potential. 

This reduces the three-dimensional NLSE (|l|) to: 



^pf{x)h{y,z) (3) 

Eq. (H) may be projected onto the ground state of 
h[y, z) and integrated over the transverse dimensions y, z: 



^^^^(y'^)[^-2M^ + LAt 

+V{x,y,z)]f{x)h{y,z)^0 (4) 



As will be shown in Sec. IIIB, in the limit that 



Ly, L 



J ^, z) takes the form of the ground 
state linear quantum mechanics particle-in-a-box solu- 
tion hgs{y,z) = hoSin{Try/Ly)sin{'!Tz/Lz). Requiring 
hgs{y,z) to be normalized to 1, ho = 2. We multiply 
Eq. ^ through by 2M^^/?i^ and obtain a simple, quasi- 
lD NLSE: 
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^V{x,y,z)\j{x)^^ 



(5) 



where / is the dimensionless wavefunction describing ex- 
citations along L; | / P /L is the longitudinal part of 
the single particle density; V(x) = {2M^ ITi')V {x) is 
the confining potential; and p = {2M^'^ /h?)p is a dimen- 
sionless chemical potential which is now the eigenvalue 
of the problem. 

The notation is further simplified by combining the 
longitudinal length of the confining potential and the 
healing length into a single dimensionless scaling param- 
eter: 



(6) 



A is an important parameter which will determine many 
of the properties of the stationary states. For the purpose 
of mathematical ease, x = x/L. Using the approxima- 
tions for Ly , Lz and dividing through by the integration 
factor of 9/4 results in the dimensionless, ID NLSE we 
shall use henceforth: 



[~>?effdl+\f{i)? +V{i)]f{i) = Peff f{i) 



(7) 

where jl^ff — fi — Stt-^/Q and A^^^ = 4A^/9. For the 
purposes of this presentation = 4/(9 * Stt \a \ p) and 
p = {2M(^'^ /h^)p — Stt^/Q. However, we shall simply 
drop the eff subscripts, as such constant factors make no 
difference in our results. 

For comparison with experiment the conversion factors 
from the dimensionless p to p in pK are listed. The gen- 
eral conversion is ^ = (8.34 * 10^^^){pa/M)p,, where M 



IS m atomic mass units, p is m cr 



and a is in nm. Us- 



ing common experimental values^ oi p ^ 10^'*, for ^^Na 
a ~ 2.75, and for ^^Rb a ~ 5.77, the conversion factors 
are 0.0723 and 0.0401, respectively. Since the dimension- 
less chemical potentials found will be on the order of 1 to 
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10 this gives a sense of the energy scale of the solutions, 
on the order of 0.1 to 1 fj,K. Note that throughout the 
presentation an experimentally reasonable test scale of 
A = 1/25 will be used for illustrative purposes. 

As I /(i) p is a single particle density, it is normalized 
to 1 rather than N: 



dx\fix)\^^l 



(8) 



The number of atoms N, which is proportional to the 
coefficient to the nonlinear term in Eq. (|^), is then con- 
tained in the ratio of the healing length to the box length, 
AociV-i/2. 

The NLSE (0), subject to normalization (^) and such 
boundary conditions as will be described below, is the 
equation we will solve. 



III. BOX BOUNDARY CONDITIONS 

We now consider the solution of Eq. (|^) in regions of 
constant potential, which may be taken to be V{x) = 
without loss of generality. We note first that if f{x) van- 
ishes anywhere in an interval, as for example at the edges 
of the box, then /(i) may be taken to be purely real 
throughout that interval. This is easily established by 
considering a Taylor series expansion of / in the neigh- 
borhood of the point at which it vanishes. Thus, we may 
remove the absolute value symbol in Eq. (0) and so re- 
cover an ordinary nonlinear equation for a real function: 



(9) 



By multiplying through by /, integrating both sides, 
and then solving for dx, the solution may be written in 
the form: 



V2A 



f/jR- 



dt 



(10) 



where R± = IzbVl ~ C, with C a constant of inte grat ion, 
and m = R-/R+. Comparing Eq. ( |lo|) to Eq. (Al) in 
App. it is apparent that they differ only by trivial 
scaling factors. Therefore in the box the most general 
solution is a Jacobian elliptic function, which as shown 
in App. ^ must be the sn function. A brief review of the 
form and properties of the Jacobian elliptic functions is 
given in the appendix. 



A. Solutions and spectra 

The most general form of the solution is: 

f{i)=Asn{kS: + 6\m) (11) 

where the notation sn(a; | m) is standard, as used in 
App. ^. k and S will be determined by the boundary 



conditions below while A and m will be determined by 
substitution of Eq. (^l]) into the NLSE and by normal- 
ization. 

The boundary conditions are: 



/(O) = /(I) = 



(12) 



The boundary condition at the origin can be satisfied 
most easily by taking S = 0. The function sn(x | m) is 
periodic in x, with period equal to 4i^(m), where K{m) 
is an elliptic integral of the first kind (see App. A). Thus 
the boundary equations at a; = and i = 1 are satisfied 
if fc = 2jK{m), where j G {1,2,3,...}. The number 
of nodes in the jth solution is j — 1. We will give a 
more general interpretation of j below. We then solve 
Eq. (0) by substituting Eq. (|ll|), using Jacobian elliptic 
identities, and setting coefhcients of equal powers of sn 
equal. This results in equations for the amplitude. A, 
and the chemical potential, fl: 



A^ = 2m{2jK{rn)f\^ 



(13) 
(14) 



Substituting Eq. ( p^ ) into Eq. (^), and noting that the 
integral over sn^ can be defined in multiples of the quarter 
period K{m), we obtain the normalization condition: 



2{2jK{m)fX\l-§P\)^l 
Kim) 



(15) 



where E{m) is the complete elliptic integral of the second 
kind. Eq. (O) then becomes: 



f{x) = V2m{2iK{m))Xsi\{2iK(m)x \ m) 



(16) 



This leaves the chemical potential (|lj) and the wave- 
function ( p^ determined up to the parameter m and 
the scale A. In Fig. |l| a graphical solution of Eq. (|l5| ) 
is shown. The plot demonstrates that the solutions are 
unique. Such solutions are in one-to-one correspondence 
with those of the ID particle-in-a-box problem in linear 
quantum mechanics. 

Plots of the wavefunction for the ground state and the 
first three excited states are shown in Fig. ||. To meet 
the box boundary conditions the wavefunction drops to 
zero over the scale of the healing length ^. When the 
zeros of / are well-separated, the analytic behaviour of / 
near a zero, Xo, approaches / ~ tanh((i — io)/(A\/2)) = 
ta,nh{{x — Xo)/ {(,V2)). We refer to this behaviour at each 
node of / as a kink. is proportional to the density of 
particles in a BEC; this density is constant everywhere 
except at the boundaries and the kinks, where it dips to 
zero. 

In Fig. H we plot the chemical potential spectrum of 
this solution type as a function of A~^, the number of 
healing lengths per box length. The leftmost portion 
of the plot corresponds to the particle-in-a-box limit and 
the rightmost portion to the topological soliton limit. We 
now discuss these two limits. 
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B. Particle-in-a-box limit 

High chemical potential states in which the kinks over- 
lap become particle-in-a-box type solutions, as can be 
seen in Fig. gd. This is both the zero density, linear 
limit and the highly-excited-state limit. Mathematically, 
m ~f 0^ and sn — > sin. Physically, j A ^ 1. In this limit 
K{m) TTn/2 + m/8 + 0{m^)) and m l/O'TrA)^, so 
that Eq. ( [14| ) becomes: 



/-i-jVA2(l + 



3 1 



(17) 



corresponding to: 



M = 



J'TT 



(1 



3to 



2_2t2 



J'TT 



2ML' 



-(1 



12aNL 
Atfn 



(18) 



which clearly converges to the well-known linear quantum 
mechanics particle-in-a-box chemical potential as m — > 
0+. 

One may also obtain this result from first order per- 
turbation theory. The Hamiltonian for Eq. (|^) is = 
Ho + Hi, where Ho = ~{h^/2M)dl Hi = {g/{LAt))P, 
and the box boundary conditions are implicit. As the 
solutions are real, we have dropped the absolute value 
sign in Hi. Note that we have put the units back in. 
The solution to Hof = A^^*^'/ is / = sin(7rji) with 
/i(o) = {fn^h'^)/{2ML'^). The first order perturbation- 
correction to /i^*^^ yields: 



(1) 



AnhaN 1 3 
M 



(19) 



where we have substituted in the definition of the cou- 
pling constant g. By noting that = 2 in Eq. (|l3 ) 
in this limit, and using Eq. ( p^ ) together with Eq. ^) 
to eliminate Na in favor of the parameter m, one re- 
covers the same first order perturbation-correction as in 
Eq. (F 



^ ~ 2ML^ T" 



(20) 



C. Topological soliton limit 

One may add a kink without disturbing another kink, 
provided that the overlap between them is exponentially 
small in the ratio of their separation to the healing length. 
In analogy with vortices the chemical potentials of the 
kinks ought to be additive. This is thcj-.large particle 
number, highly nonlinear, Thomas-FermicI limit. Math- 
ematically, m — > 1^ and sn tanh, formally called a 



topological soliton. Physically, N ^ 00 implies A ^ O"*", 
so that we expect the kinks should never overlap. We 
note that in a BEC experiment the box length is held 
fixed while atoms condense. 

By solving for K{m) in Eq. ( p^ ) and using E{m) 1+ 
(see App. |A|) we find that K{m) — )■ k where: 



so that Eq. ( |l^ becomes: 

M = 2{2jfX'^K^ 
while Eq. ( p^ ) becomes: 

f{x) = V2{2jK)Xsii(2jKx I to) 



(21) 



(22) 



(23) 



We have found that this limit suffices to calculate 
chemical potentials for which j < (1/(5A)) to better than 
1%. This estimate assumes an overall scale size of ~ 5f 
per kink. The chemical potentials for the j = 1, j = 2, 
and J = 3 solutions shown in Fig. ^ satisfy this crite- 
rion, for example, as do any ground states for a healing 
length of smaller than 1/10. If we now further require 
that (SfX^) < 1 then Eq. (|2|) becomes: 



/2 = 1 + 2V2jX 



(24) 



from which we see that the chemical potentials of addi- 
tional kinks are indeed additive. This additivity is ap- 
parent in Fig. I in the Hmit A^^ = L/^ 00. Eq. (|2|) 
is identical in form to that of the harmonic oscillator in 
linear quantum mechanics. 

Putting back in the units of h^/{2M^'^), we find that 
the chemical potential for formation of an additional 
kink, A/Lt, is proportional to ^/N : 



(25) 



so that the chemical potential to add a kink increases as 
atoms condense, when all units are included. 

We may also solve for the excitation energy to add 
an isolated kink to the N body system by findijig the 
expectation value of the many-body HamiltonianH: 



(H) 



iV[A2 



1/ 



' |2 



(26) 



where just as for the chemical potential, (H) = 
{2M / 1)?) {H) is a dimcnsionless energy. Substituting 
in the stationary solution Eq. (|ll|), one finds: 



{H) = N—j^KimfX^l (1 + 2m)K{m) - (1 + m)E{m) ] 
3 

(27) 

which in the limit as (8j^A^) <C 1 becomes: 



4 



(28) 



Therefore the excitation energy to add a kink is 
A{H) ~ 2N^/2X, or putting the units back in: 



2o3 



h^2 



2M 



naAt _3/2 



(29) 



This is just as one would expect; the energy to add a kink 
no longer depends on the box length when the kink size 
is much smaller than the box. Instead it simply depends 
on the density. 



IV. PERIODIC BOUNDARY CONDITIONS 

There are three solution types for periodic bound- 
ary conditions. There are constant amplitude solutions 
which are plane waves; real symmetry-breaking solutions 
similiar to those found in Sec. 
complex symmetry-breaking so. 



Ill; and a novel class of 
utions. The former two 



are in one-to-one correspondence with particle-on-a-ring 
solutions in linear quantum mechanics, while the latter 
one is only found in the presence of nonlinearity. As the 
ring is rotationally invariant, the symmetry-breaking so- 
lutions will have a high-degeneracy, in analogy with vor- 
tices in two dimensionsEa. The periodic boundary condi- 
tions are: 

/(O) = /(I) (30) 
/'(O) = /'(I) (31) 



A. Constant amplitude solutions 

If we assume that the amplitude is constant then we 
obtain plane wave solutions of the form: 



(32) 



where n £ {0,±1,±2, ...}. The amplitude is constrained 
by normalization to be 1. Substituting Eq. ( |3^ ) into 
Eq. (^ we find the chemical potential: 



A27m f 



(33) 



from which we obtain the lower limit of the chemical po- 
tentials under periodic boundary conditions, /i = 1. This 
is just what we expect physically for the repulsive BEC. 
The ground state on a ring is the condensate spread out 
evenly. There is no symmetry breaking. For n 7^ each 
solution is two-fold degenerate, as n can be either posi- 
tive or negative, while the n = 0, ground state solution 
is non-degenerate. 

Note that these states could also be termed angular 
momentum eigenstates or quantized rvortices, as for ex- 
ample in the work of Matthews et alB3. 



B. Real symmetry- breaking solutions 

As we have exchanged the ring for the box, Eq. (|ll]) is 
the real solution. One simply changes k from 2jK{m) to 
4:jK{'m) in order to satisfy Eqs.(^0|) and (^l|), i.e. from 
multiples of the half period to multiples of the whole 
period. The number of nodes will be 2j rather than j-1, 
where j e {1,2,3, ...}. We temporarily keep S set to 0. 
But note that, unlike for box boundary conditions, under 
periodic boundary conditions S is arbitrary. 

Then all the results from section III hold with the new 
fc, by letting j — > 2j in all equations. The energy and 
wavefunction are determined uniquely by graphical so- 
lution of Fig. |l|. In Fig. || we show the first two states. 
Both the linear quantum mechanics, particle-on-a-ring 
limit and the topological soliton limits are reproduced. 
In the latter the same kind of non-overlapping criterion 
applies as before. Thus given (16j^A^) ^ 1 Eq. (|lj) 
becomes: 



= 1 + 2\/2(2j)A 



(34) 



Note that the factor of two in front of j shows that, on 
a ring, kinks of this type come in pairs. 

If S is permitted to vary arbitrarily, the degeneracy in- 
herent in these symmetry-breaking solutions is obtained. 
The entropy associated with j kinks depends logarithmi- 
cally on the box length L, and, since there are approx- 
imately possible positions for the kink, the entropy 
is: 



S kBln{ 



1 



(35) 



where the factor of 4-^2 comes from 2V2C for each of the 
two kinks. This is consistent with the non-overlapping 
criterion we used in obtaining Eq. (E4). 



C. Complex symmetry- breaking solutions 

For complex solutions we divide the wavefunction into 
a phase and amplitude: 



(36) 



and obtain four boundary conditions. From substituting 
Eq. (|6|) into Eq. ( pO[ ) and taking real and imaginary 
parts: 



r(0) = r(l) 



(1) - 0(0) = 27™ 



(37) 



(38) 



where n is an integer which we will call the phase quan- 
tum number. From substituting Eq. ( ^6|) into Eq. ( [3l| ) 
and again taking real and imaginary parts: 



r cos( 



rm smc 



\x=o- 



\r cos( 



sm ( 



U=i (39) 



5 



[r' sin (j) — rcj)' cos cj)] \x=q= [r' sin cf) — rcf)' cos 0] \x=i (40) 

Substituting Eq. (|6|) into Eq. (@), we divide the NLSE 
into real and imaginary parts. We integrate once to solve 
for 0' in the imaginary part, and find: 



a 
S 



(41) 



where a is an undetermined constant of integration and 
S = r{xY is the single particle density |/(x)p. Substi- 
tuting Eq. ( pi] ) into the real part, we multiply through 
by r' and integrate again. We find: 



(5') 



l\2 



-2[- 



1 

A2 



2M c.2 



A2 



(42) 



where (3 is an additional undetermined constant of in- 
tegration. A similiar solution method has been used 
by Drazin and Johnson in an elementary discussion of 
solitonsEl. For the complex solutions, Eqs. (p])-(^2[) re- 
place the NLSE as the equations to solve, together with 
boundary conditions (p7|)-(p0|), and the normalization 



We may rewrite Eq. (B2h as an integral: 



dS 



2v/(A-2)53 + (-2/iA-2)S2 + pS- 2a2 



(43) 



This is an elliptic integral. Any elliptic integral can 
be expressed as the sum of a finite number of elliptic 
integrals of the first, second, and third kinds. Given that 
a and [3 are real, such integrals may be reduced to a 
standard form with Cayley transformations, so that < 
m < 1 and the parameter m is realEj. Therefore all 
intrinsically complex solutions to the ID NLSE may be 
written as a sum over standard elliptic integrals. 

We have found real symmetry-breaking solutions for 
which the density is proportional to sn^ . These solutions 
vanish at 2j points around the ring. We look for solu- 
tions of a similiar form for which the density does not 
vanish. The physical motivation for such a solution type 
will become clear in Sec. Using our physical intuition, 
we are able to bypass the use of Cayley tr ansf ormations. 

From the Jacobian elliptic identity Eq. ( A2): 



sv?{x I m) = — (1 — ch\^{x I m)) 
m 



(44) 



where dn(5; | m) is the Jacobian elliptic function which 
we describe in Fig. 0. We thus generalize the real 
symmetry-breaking solutions, Eq. ([L1[), as follows: 



r2(i) = A^{1 + jdn'^[kx + S \ to)) 



(45) 



where — 1 < 7 < 0. 

By setting k equal to the full period of dn, i.e. k = 
2jK{m), and S = 0, we will automatically match the 
boundary conditions related to amplitude, Eqs. (37), 
(p9|) , and (EO) . It will remain to satisfy the phase quanti- 



to the number of nodes or kinks. Here, as the density no 
longer goes to zero, j is to be interpreted as the number 
of dips, or density- notches as we will call them, in the 
density r(i)2. This is consistent with our previous defi- 
nitions of j. 7 is then the depth of the notch, while 
is put in to satisfy normalization. We will consider the 
case of general 6, and thus degeneracy, later. 
Eq. (Ea) then becomes: 



r^{x) = A-^ (l+jdn-^ {2 jK{m)x \ to)) 



(46) 



Substituting this into Eq. (^), using additional Jaco- 
bian elliptic identities, and setting coefficients of equal 
powers of dn equal, we obtain four equations in the pa- 
rameters a, /3, 7, and k. Eliminating we are left with 
a, 7, and A^ as a function of m. A, j, and p,. We use the 
normalization (^) to constrain jl and find: 

/i = ^ + I2f\^ E{ni)K{m) - Af{2 - m)\^K{mf (47) 

From this we obtain our equations for the parameters a, 
7, and A^: 

-8/(1 - m)\^K{mf) * (1 - %f\^K{mf 
+Mf\'^E{mfK{mf - lQpX^E{m)K{m) 
*{-l + A3^\^K{mf)))f/^ (48) 



1 + %pX^K{m)E{m) 



l + 8fX^K{Tn)E{m) 



(49) 



(50) 



zation (38). For real solutions we said that j was related 



This leaves the constant of integration a in cj)' — a/ r^, 
the depth 7, the prefactor to the density A^, and the 
chemical potential /i, determined up to the number of 
density-notches j, the scale A, and the parameter to. For 
a given A and j we then numerically integrate the phase 
(|l|) and, using TO as our free parameter, adjust to until 
the boundary condition (38) is met, i.e. until the phase 
quantum number n is an integer. We note that all pa- 
rameters are monotonic in to so that our algorithm is 
quite straightforward. By symmetry of the ring, n can 
be either positive or negative, so that each solution is 
two-fold degenerate, just as we found for the constant 
amplitude solutions. 

In Fig. U we show the amplitude and phase of one and 
two density-notch solutions at our test scale of A = 1/25. 
We have plotted the amplitude above the phase to make 
apparent that the phase is a background constant slope 
with a region of increased slope where the density-notch 
occurs. The deeper the notch, the larger the increase in 
slope. In the limit that the notch dips to zero to form a 
node, the phase becomes a step function of height tt per 
step and the real solutions are recovered. 
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If S is generalized so that it is arbitrary, a similiar 
degeneracy to what was found in Eq. (|35| ) results: 



S ^ kBln{——) 
OA 



(51) 



where C = 7r\/6, as we shall show in Sec. [V D 



D. Bounds 

Real stationary states can have an arbitrarily large 
number of nodes. But the number of notches for nodeless 
states is limited. We set three bounds on the complex, 
nodeless solutions; the maximum chemical potential; the 
minimum and maximum phase quantum number; and 
the minimum scale to obtain j notches. As a conse- 
quence of these bounds there are some scales at which 
no complex solutions exist. 

The maximum number of density-notches that can fit 
on the ring is obtained from the lower limit on the period 
of the dn function in Eq. (Q). When the notches over- 
lap too much they are no longer solutions to the NLSE. 
The dn function approaches its minimum period of tt as 
TO — > 0+ . In this limit Eq. (j^j is the maximum chemical 
potential: 



(52) 



In this same limit the amplitude approaches a con- 
stant which the normalization constrains to be 1. From 
Eqs. (|4^) and ( ^8| ) we find a relation between the max- 
imum number of density-notches, the phase quantum 
number, and the scale A: 



A = 



1 



(53) 



One may invert this relation to find the maximum n 
for a given number of density-notches j. Eq. ( p3| ) requires 
that n > j /2. Thus there is both an upper and a lower 
bound on n: 



J 



1 , 1 



8V2A2 



(54) 



As j increases rimax ^ {j f'^)^ ■ Only integer n can solve 
the phase quantization condition (p8|). It follows that, for 
a given A"-'^, more odd j solutions will be available than 
even j solutions, because j/2 for odd j is half integer. 
This is apparent in Fig. ^ 

The above bounds imply that at a given scale the num- 
ber of density notches is bounded from above. In Fig. ^ 
we plot the scale at which j density-notches become pos- 
sible. For less than 7.7 healing lengths to the box length, 
there are no complex solutions at all. Then, in order, 
j = 1, 3, 5, 2, 7, 9, 4, ... solutions become possible: 



a;^ 



mm 
min 



7rV2 + A jo. 



dd 



(55) 



for even and odd j, respectively 

The minimum inverse scale for j = 1 is A^^ — 7r-\/6. 
This is the natural size of a density-notch. At smaller 
inverse scales the notches are affected by overlap, so 
that complex solutions do not exist, while real solutions 
become sinusoidal. Thus the parameter C in Eq. ( [5l| ) 
takes the value C, ~ tta/G. Checking the limits, as 
A~^ {\~^)min — Ti'VS, 5^0+. There is only one 
configuration. 



E. Spectra 

We show the chemical potential spectra as a function 
of A^^ for the three types of stationary states on the 
ring: real, constant amplitude, and intrinsically complex. 
In Fig. ^ the two lowest real spectra are shown. For 
comparison we have overlayed the four lowest constant 
amplitude spectra on the same figure. In Fig. |^ we show 
the three lowest spectra for the complex solutions. 

For our experimentally reasonable test scale of A = 
1/25 the order of stationary states, starting with the 
ground state, is: constant amplitude, singly quantized 
vortex, single grey density notch, doubly quantized vor- 
tex, real two-node solution, two grey density notch, etc. 
Note that the minimum chemical potential la jl — \. 

Since the real solutions are a limiting case of the com- 
plex solutions, the two scale in the same way and their 
energy levels do not cross. But the constant amplitude 
solutions depend differently on inverse scale, so their en- 
ergy levels can cross with those of the other solutions. 



V. CONNECTION WITH SOLITON THEORY 

The dimensionless, time-dependent, free NLSE is: 

[ludt + ed.,:- I fix, t) p]/(a:, t) = (56) 

where v = jU has units of time and we have chosen 

to use X rather than x. The single grey or dark, density- 
notch solution to this equation takes the formed: 

fix-ct,t) = V2[t^ + S,Ktanh{K{x - ct))]e-'''^/'' 



2^ 



(57) 



where k is the width, c is the speed, and ^ is the chemical 
potential. 

Using an identity to change from tanh to sech, the 
amplitude squared is: 

\f{x -ct) p= 7/ - 2^^K^sech^{K{x - ct)) (58) 

where 77 = {{c^iy^)/ {2S_'^) + 2k^(^). The speed is: 
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(59) 



where c varies from zero to the Bogohubov sound speed0 

which, on insertion of constants, is Cmax — \J A'Kh'^ ap / m? . 
For c = we recover the dark, tanh soHton mentioned in 
Sec. [II. For c c^ax the grey density notch approaches 
zero depth. Thus the speed of the density notch has 
a finite range and the zero depth solitons approach the 
speed of sound from below. 

The dn function is a periodic generahzation of a sech. 
Therefore replacing the sech in Eq. ( |57| ) by a dn is equiv- 
alent to replacing a soliton by a soliton train. With 
this substitution, Eq. ( [5^ ) becomes identical in form to 
Eq. p^ ) up to its time dependence. A single soliton is set 
in motion by a phase jumplHlo. The same holds true for 
a soliton train. In the plots of the phase of the complex 
solutions. Figs. ^ and there was a background linear 
slope with a jump where the density-notches occurred. 
Thus the complex stationary states can be interpreted as 
momentum boosts of the condensate with grey soliton- 
trains superimposed. The velocity of the momentum 
boost exactly cancels the velocity of the soliton-train, 
resulting in a stationary state in the lab frame. 

Although free grey soliton-trains can vary continuously 
in speed up to Cmaxj periodic boundary conditions require 
that the boost speed which brings about our stationary 
solutions is quantized. Using Eq. ( |59| ) and making the 
substitutions rj — A and —2k^ = Aj so that we change 
over to the constants used in Eq. (Eq), the boost speed 
of the condensate which results in cancellation is: 



Cboost = -v/2A(l+7) 

V 



(60) 



where 7 — > —1+ is the c = limit and 7 = is the 
Bogoluibov sound speed. 7^—1+ repro duces the real, 
dark soliton solutions found in Sec. IV B . Thus the real 
solutions are soliton-trains at rest. 



VI. CONCLUSION 

We have presented the complete set of stationary so- 
lutions to the nonlinear Schrodinger equation under pe- 
riodic and box boundary conditions in one dimension for 
the case of repulsive nonlinearity. In a box all solutions 
may be taken to be real. On a ring there are three solu- 
tion types: constant amplitude solutions which are plane 
waves; real symmetry-breaking solutions; and a novel 
class of complex symmetry-breaking solutions which cor- 
respond to a boost of the condensate in one direction 
with density-notches moving in the opposite direction, 
so that they are stationary in the lab frame. 

Real and constant amplitude solutions are in one-to- 
one correspondence with those of the analogous particle- 
in-a-box and particle-on-a-ring problems in linear quan- 
tum mechanics. Complex, symmetry-breaking solutions 



are uniquely nonlinear. We showed that solutions of 
non-constant amplitude may be treated as density-notch 
soliton-trains. As the natural size of a density- notch is 
7rV6, the minimum scale size needed to obtain complex 
solutions is i/i^ — tt^/G. In the context of the BEC, this 
means that the number of atoms determines the available 
solution-types. 

Th&|p3sults of this paper cast the findings of previous 



workMcd in the larger context of a comprehensive set of 
solutions to the NLSE. In the previous work, perturba- 
tions of the phase of the purely real, box-type solutions 
were found to induce solitonic motion. Such perturba- 
tions have recently; been applied in the laboratory to, 
three-dimensionalli3 and nearly quasi-one-dimensionalll3 
BECs and have generated solitary waves. Our present 
work gives a larger class of solutions, whose response to 
perturbations may suggest further experiments. 
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APPENDIX A: JACOBIAN ELLIPTIC 
FUNCTIONS 

We here briefly review the properties of Jacobian el- 
liptic functions and establish the notation used herein. 
There are a total of twelve such functions. All of them 
solve the unbounded ID NLSE in one form or another. 
Of the twelve, six are normalizable. Of these, only three 
have different physical form. They are the sn, the cn, 
and the dn. We plot them in Fig. ^. Of these, only the 
sn solves the NLSE for repulsive nonlinearity while only 
the cn and dn solve the NLSE for attractive nonlinearity. 
The cd also solves the repulsive case. It differs from the 
sn by a shift of a quarter period but is otherwise identical. 

The function sn(u | m) may be written in integral form: 



dt 



VT^VT 



(Al) 



where u = sn~^(x) so that x — sn{u \ m). The functions 
cn(u I m) and dn(u | m) may then be defined by the 
equations: 



cn(u I m) = ^ 1 — dn^ [u \ m) 
dn(it I m) = — msn^(M | m) 



(A2) 



The limits of the sn, cn, and dn functions, along with 
the complete elliptic integrals K(m) and E{m) are listed 
in Table ^ The period of the sn and cn functions is 
4K{m), while that of the dn function is 2K(rn). We 
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directpthe reader to Abramowitz and StegunE^I and other 
worksEj for a further review of the properties of these 
functions. 
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FIG. 1. This graphical solution of Eq. ( |15[ ) shows that for 
a given scale and number of nodes the real solution to the 
stationary NLSE under box or periodic boundary conditions 
is unique. A is the scale and j — 1 with j € {1, 2, 3, ...} or j 
with j G {2, 4, 6, ...} is the number of nodes, respectively. The 
three curved lines are plots of Eq. (^ solved for the number 
of nodes j, with A"^ = L/C = 10, 25, 50. The left-hand side of 
the plot is the m — 0, linear limit, while the right-hand side 
exponentially approaches the m = 1, topological soliton limit. 
The solutions are found where these lines intersect with the 
horizontal lines of j. Note the rapid convergence to m = in 
the high j limit, so that for large j the solutions are in the 
linear regime. 
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FIG. 2. Real stationary solutions to the NLSE under box 
and periodic boundary conditions. These are in one-to-one 
correspondence with those of the analogous particle-in-a-box 
and particle-on-a-ring problems in linear quantum mechan- 
ics, and may also be characterized as dark soliton trains. 
Box: (a)-(d) are the ground state and first three excited 
states. Ring: (b) and (d) are the first two solutions of this 
type. Chemical potentials: (a) fl — 1.120 (b) jl — 1.253 (c) 
fi = 1.402 (d) /i = 3.028. All plots are for the test scale of 
= 1/25. 



FIG. 4. Complex, symmetry-breaking stationary solutions 
of the NLSE on a ring. These are grey density-notch soliton 
trains. A supercurrent around the ring exactly cancels their 
velocities to make them stationary states in the lab frame. 
j is the number of density-notches, n is the phase quantum 
number, and jl is the chemical potential, (a) Amplitude and 
(b) phase/27r of the j = 1, n = 1, p. = 1.197 stationary 
state, (c) Amplitude and (d) phase/27r of the j = 2, n = 2, 
p, — 1.331 stationary state. All plots are for the test scale of 
^/L = 1/25. 
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FIG. 3. Chemical potential spectra of real stationary 
states, as a function of inverse scale L/^, with stationary plane 
wave spectra shown for comparison. Solid lines: shown are 
n = 0, 1,2,3, where n is the phase quantum-number of the 
plane wave on the ring. Dashed lines: real stationary states 
of the NLSE in a box and on a ring are soliton-trains. Shown 
are j — 1, 2, 3, 4 with j — 1 the number of nodes in a box and 
j — 2,4: the number of nodes on a ring. Note that for very 
fine scale, i.e. L/^ large, the chemical potentials are evenly 
spaced. This corresponds to the topological soliton limit in 
which the chemical potentials are additive, just as for vortices. 




FIG. 5. Chemical potential spectra for complex, symme- 
try-breaking stationary states on the ring as a function of 
inverse scale. The three solutions shown here we re fou nd us- 
ing the numerical algorithm prescribed in section IV C . They 
are {j, n) — (1, 1), (2, 2), (3, 2), where j is the number of den- 
sity-notches and n is the phase quantum number. The left- 
most points are the minimum inverse scale and maximum 
chemical potential possible for such a solution. Bounds and 
ordering of the solutions are shown in Fig. pi 
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FIG. 6. Minimum inverse scale for j density-notches to 
become available. The lower curve is odd j and the upper 
curve is even j. Note that at a given inverse scale there are 
many more odd solutions than even solutions available. The 
ordering of the solutions is j = (1, 3, 5, 2, 7, 9, 4, ...). 
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FIG. 7. The three basic Jacobian elliptic functions, for 
the parameter m = 0.99. The solid line is sn, the dotted 
line is cn, and the dashed line is dn. Note that the period 
of dn is half that of the other two and that dn is nodeless. 
All Jacobian elliptic functions may be constructed from these 
three. Of the twelve possible functions, these three shapes are 
the only normalizable ones not different from each other by 
a translation along the horizontal axis or a renormalization 
along the vertical axis. 



TABLSi I. Limits of Jacobian elliptic functions and 
integralsCJ. 



m = m = 1 

sn(u I m) sin(it) tanh(ii) 

cn(w I m) cos(ii) sech('u) 

dn(-u I m) 1 sech(tt) 

K(m) 7r/2 oo 

E{m) Ti/2 1 



11 



